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The ratio of the muon flux density to charged particle flux density at distances of 300 and 600 m 
' from the shower axis (p^ (300) / ps (300) and p^ (600) / ps (600)) is measured. In addition, the energy 

dependence of p^ (1000) is analysed for showers with energies above 10^* eV. A comparison between 
the expeiimental data and calculations performed with the QGSJET model is given for the cases of 
^ ^1 piimary proton, iron nucleus and gamma- ray. We conclude that the showers with Eo > 3 X 10^* eV 

I ' can be formed by light nuclei with a pronounced fraction of protons and helium nuclei. It is not 

O ' excluded however that a small part of showers with energies above 10^® eV could be initiated by 

primary gamma-rays. 

■ Keywords: Cosmic Rays; Extensive Air Showers; Mass Composition. 



o 



X 
J3 



The Yakutsk Complex EAS array is able to record muons with the threshold energy 

£;th IGeV X sec(6i). 

For this purpose the array is equipped with 6 muon detectors with the full area of about 
300 rv?. The detectors are located at the distances R ^ 100, 300, 500 and 1000 m from the 
array center ^ The muon detectors are interrogated upon the availability of the array "mas- 
ter", i.e. the coincidence of responses from the scintillation detectors of three observation 
stations located at the top of 500 m arm equilateral triangles. After almost 20 years of con- 
tinuous work of the array, a high statistics on the showers with energies above 10^^ eV has 
been accumulated. 

For the present study, the Yakutsk array data accumulated during the period from 1994 
to 2004 (i.e. after the array modernization^) was taken into consideration. Thus the ac- 
curacy of the EAS axis determination has been increased twice. This considerably im- 
proved the quality of the analyzed data. The showers with £^0 > 10^^ eV and zenith angles 
< 60° have been chosen for the analysis. It was required for the shower axis to be within 
the boundaries of the array and that at least three muon detectors (one of which is at a dis- 
tance of 1000 m from the axis) response to the shower The muon flux density, p^ (lOOO), 
is then determined as the median between the densities adjusted to (B) = 1000 m in accor- 
dance to the mean muon lateral distribution function. For this purpose, besides the detectors 
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Fig. 1 . a) Fraction of muons with energies above 1 GeV (%): (300) /ps (300) (•) and (600) / ps (600) (o). 
b) Pp (1000) as a function of Eq for tlie observed events at energies witliin tlie ranges 10^® - lO^'' eV (squares), 
10^^ - 10^" eV (filled circles) and 10^° - 10^^ eV (triangles). The expected ±1<t limits for the distributions 
are indicated (in both panels) by different curves for the proton, iron and gamma-ray primaries, as is explained in 
the text. 
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responded in the shower, the data from the threshold detectors were also included into the 
analysis. The muon flux densities at the distances 300 and 600 m from the shower axis 
were calculated similar way. The results of our study are shown in Figs. 1 a and b. 

Calculations suggest^ ^ that the function [R) slowly depends on the zenith angle in 
the angular interval /S.9 — (0° — 60°). Hence, within this interval it will be is a reasonable 
approximation to apply the mean angle 9 = 39° for a comparison of calculations and 
experimental data. In order to take into account both physical and methodical fluctuations 
in the measurements of p^(lOOO) (including the implementation errors and the errors in 
determination of the shower axis location) we considered the "judicial" interval in one 

—mean ^phys _i_ _meth 

In Figs. 1 a and 1 b the results of calculations are shown by dotted lines for the case 
of primary proton and by dashed lines for iron nucleus. Solid lines indicate the upper and 
lower limits for the EAS with Eq > 10^® eV initiated by a primary 7-quantum. In this case 
we adopted the results of calculations from Ref.'T' The dots and circles in Fig. 1 a are for 
the experimental data on (300) and (600), respectively. In Fig. 1 b, the showers in the 
energy range of 10^^ — 10^^ eV are shown by squares, the showers with Eo > 10^^ eV are 
shown by dots, and the showers of maximum energy are shown by triangles. 

Our comparison of experimental data presented in Fig. 1 a with the calculations car- 
ried out for the case of the primary proton and iron nucleus and based upon the QGSJET 
model, confirms the hypothesis that a considerable fraction of ultimate energy EAS is be- 
ing formed by protons. The proton fraction of the primary cosmic rays decreases below 
the energy of about 10^^ eV. It can be seen from Fig. 1 b that the basic amount of points 
relevant to the showers having energies larger than 10^^ eV falls into the "proton inter- 
val"; 23 points fall into the zone of superposition of primary protons and iron nuclei; and 
only 19 showers among 116 fall into the zone near the upper boundary evaluated under 
assumption that the primary particle is a 7-quantum. Our analysis indicates that, even if 
the fluctuations were taken into account, a substantial probability remains that the showers 
of such energies are generated by neutral particles, in particular, by primary 7- quanta. The 
analysis of arrival directions of showers with Eq > 10^^ eV is justified for searching the 
sources of highest energy cosmic rays. The accuracy of determination of the arrival angles 
for these showers is expected to be not worse than (0.5° — 1.5°). 
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